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Abstract 

We describe some properties of nuclear ground states using a microscopic Vlasov 
approach. We propose a new method to find the energy minimum of the system 
enforcing the fermionic nature of the system. Calculations are performed for two dif- 
ferent parametrization of the mean field: a simple momentum independent Skyrme 
potential and a momentum dependent one as proposed by Gale, Bertsch and Das 
Gupta. We show that the binding energies and radii of nuclei from ^^O to ^''^Pb 
are reasonably reproduced by opportunely tuning the surface term. We also made 
some calculations for neutron rich/poor nuclei. 



1 Introduction 



Microscopic Vlasov simulations have been applied in the last two decades or 
so with varying degree of success to heavy ion collisions (HIC)[1]. While col- 
lective effects, particles distributions and production have been successfully 
described, the approach fails when many particles correlations are important 
such as in fragmentation reactions. This is of course expected since we are deal- 
ing with a mean field approach. The method has been applied more recently to 
metallic cluster collisions and quark dynamics as well [2]. One of the reasons 
of the success of the Vlasov approach is surely the simplicity of its numerical 
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solution. Furthermore the Vlasov equation (VE) might inchidc some quantum 
effects as well[3]. In that respect recall that the one body distribution function 
/(r, p, t) entering the VE can be viewed as the Wigner transform (WT) of the 
one body density matrix p(r, r', t) which fulfills the Time Dependent Hartree 
Fock (TDHF) equation. In fact it can be easily shown that the VE is the 
WT of TDHF in the limit of ^ — > 0. Thus the two approaches are practically 
equivalent at large entropies where quantum effects are not important [1]. The 
numerical simplicity is one of the advantages of the VE over TDHF the other 
is the possibility of easily including a coUision term in the VE. Such a collision 
term is of course crucial for energetic HIC. On the other hand a serious short- 
coming of the VE are the calculations of nuclear ground states. While Hartree 
Fock (HE) is a quite well developed and accepted method for calculations of 
nuclear (and other fermionic systems) ground states, there is not so far an 
accepted method to obtain the ground states self-consistently in the Vlasov 
approach. This is in part due to the semiclassical nature of the VE which 
would prefer that the systems in their ground states be solids. In practical 
calculations of the VE the test particles (tp) method, shortly described below, 
is used [1] . Initially the tp are randomly distributed in a sphere of radius R 
(the radius of the nucleus) in r-space and of radius Pp (the Fermi momen- 
tum) in p-space. Such initial conditions are then evolved in time according 
to the classical equations of motion for particles moving under the influence 
of a mean field. Because of the Fermi gas initial conditions, the system ful- 
fills initially the requirement that the occupation factor, i.e. the number of 
identical particles in a cell of size 27ih, must be less or equal to one. Such a 
condition is preserved in time because the VE fulfills the Liouville theorem [1], 
i.e. the density in phase space remains constant. While this method could be 
accepted in collisions of stable nuclei at large beam energies, it can be seriously 
doubted if applied to exotic nuclei reactions and/or at beam energies close to 
the Coulomb barrier. It is the purpose of this paper to show that the ground 
states of nuclei can be also rather well calculated using the Vlasov approach. 
This will give us the opportunity to extend the calculations to exotic nuclei 
that are a quite active field now. 

The paper is organized as follows. In section 2 we review the Vlasov equation 
and the test particles method. In section 3 we apply the method to calculate 
nuclear properties with a local Skyrme potential while section 4 is devoted to 
calculations employing a momentum dependent force. We stress the important 
role of the surface and symmetry terms to reproduce some basic features of 
nuclei. A brief summary is given in section 5. 
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2 The Vlasov Equation 



The TDHF equation of motion 

dtP=[H,p\ (1) 

in the Wigner representation can be written as a series of the form: 

9J + I • V,/ - V,C7 • V,f - (2) 

where Fq[f] is the quantum corrections term which in the hmit ^ — > is 
zero [3], which is the hmit we will consider in this work. The Ihs is the Vlasov 
term and C/(r)is the mean field, which we have written in a local form but a 
generalization to a momentum dependent one is straightforward[l]. 

For the purpose of this work we will neglect the collision term in Eq. (2) and 
note that such term will be essential when dealing with HIC. 

Numerically the VE equation (2) is solved by writing the one body distribution 
function as: 

1 ^ 

/(r,P,i) = — 5^^^(r-ri(i))^p(p-pi(i)), (3) 

where the and Qp are sharply peaked distributions (such as delta functions, 
gaussian or other simple functions), that we shall treat as delta functions. 
N — Anxp is the number of such terms, A is the mass number of the nucleus. 
Actually, N is much larger than A, so that we can say that each nucleon is rep- 
resented by ritp terms called test particles (tp). The rigorous mean field limit 
can be obtained for ritp oo where the calculations are of course numerically 
impossible, even though the numerical results converge rather quickly. Insert- 
ing Eq. (3) in the Vlasov equation gives the Hamilton equations of motions 
for the tp[l]: 



-Pi 

Pi^-^nU, (4) 

ior i — 1, - ■ ■ N. The total number of tp used in this work is more than 400000 
and the results are numerically stable. 

One important feature of TDHF is that the density matrix is given by a Slater 
determinant, i.e. the wave function must be properly antisymmetrized when 



3 



dealing with fermions[4]. This imphes the following property for the density 
matrix: 




(5) 



which is preserved in time by the TDHF equation of motion. Making the 
Wigner transform of the last equation in the same limit — > gives[4]: 



i.e. the occupation density / must be zero or one. This is a direct consequence 
of the Pauli principle which states that in a cell of size h in phase space there 
cannot be two identical particles. In our work we will distinguish protons 
and neutrons with spin up or down, thus we have four types of particles. A 
consequence of this is that, for instance for a particles the Pauli principle does 
not apply since the particles are non identical. If one distributes randomly the 
tp in phase space according to the Fermi gas model, the condition given by 
Eq. (6) is fulfilled and it will be preserved in time by the VE. But the state 
so obtained is not a state of minimum energy. Notice the analogy to the HP 
method in that one initially uses a Permi gas for the wave functions. However 
a process of minimization is then applyed which finally gives the state of 
minimum energy [4] . Such a minimization is not usually done for the VE. But 
this can be easily accomplished as follows. Pirst, for numerical reasons, we put 
a grid on the available phase space and store the occupation factor in each 
point of the grid with the initial conditions given by the Permi gas model. We 
solve a frictional equation of motion with a constraint due to the occupation 
factor for each tp (/j) as: 



where /iq is a negative frictional coefficient independent on /j, and ii{fi) is the 
/-dependent frictional coefficient. is positive when /j exceeds 1, while it 

is negative for /j < 1. This /-dependence of fi is to control the local momentum 
space. The method proposed here has been first introduced in [2] and refined 
in the Constrained Molecular Dynamics calculation of Ref . [5] . 

Por illustration we show in Pig. 1 the occupation factor averaged over all the 
tp versus time for ^^O. To show the power of the method we used initially 
a very large Permi momentum. Because of that the initial occupation factor 
is rather small but it increases quickly and stabilizes around 1 after about 



P - f, 



(6) 




'i 



(7) 



4 




200 400 
t(fm/c) 

Fig. 1. Average occupation in phase space versus time for ^^O. 

200 fm/c. Results for different nuclei and forces are quite similar to the one 
illustrated in Fig. 1. 



3 Local Mean Field 



In addition to the Coulomb interaction for protons, a popularly used parametriza- 
tion of the mean field is given by the Skyrme parametrization: 

U{p{r)) = ap(r) + bp^{r) + c,VV + c„^^V^' 

where p,Pn,p are the total, neutron and proton density respectively, po — 
0.16 fm~^. The parameters a,b and n are tuned to have a compressibility 
of 225 McV[l]. The canonical values of the parameters entering the surface 
and symmetry terms are Cg = —138 MeVfm^ and Ca = 32 MeV respectively. 
However, when using a grid in phase space a surface term is spuriously in- 
troduced. Because of that we fixed the value of Cg — 6.9 MeVfm^ by fitting 
the binding energy of ^^O, i.e. a small nucleus for which the surface is very 
important while the symmetry energy is negligible. The small value of the 
surface term coefficient should not surprise, in fact for density calculations in 
r-space we use a grid of 1 fm side. This introduces a surface term spuriously. 
Changing method of solution, for instance gaussians instead of delta function, 
will give rise to a surface term from the volume part of the mean field. Also in 
that case one must adjust the surface term to some known property of finite 
nuclei as we did here. 

In Fig. 2 we plot the binding energies versus time for four different systems. 
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Fig. 2. Binding energy per nucleon versus time for different nuclei as indicated. The 
dashed hne is the Weizsacker parametrization fitted to stable nuclei. 



The dashed hne indicates the value of the binding energy obtained from the 
Weizsacker formula [4] whose parameters are fitted to stable nuclei. We stress 
that all the terms entering the Weizsacker formula arc in principle contained 
in out approach excluding the pairing energy. We see that the BE from sys- 
tematics are reasonably reproduced by the method. 

In Fig. 3 we plot the radius versus time for the systems as above. The full hnes 
refer to neutrons and the dashed lines to protons radii respectively. A neutron 
halo is quite evident for the unstable C, while the proton radius is slightly 
larger than the neutron one for the Sn nucleus. The Pb nucleus gives similar 
radii which is due to the competing action of the symmetry and coulomb terms. 
It also seems that the neutron and protons oscillate with different frequencies. 
This is especially seen in the C case, which indicates that the protons and 
neutrons have different compressibilities. We would like to notice also that 
in the Pb case, due to the strong Coulomb force it is quite difficult to reach 
a state of minimum energy. Because of the relative small number of tp per 
nucleon (about 4000), classical fluctuations could lead the system to fission if 
the minimization procedure is followed for a longer time. 
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Fig. 3. Proton (dashed lines) and neutron (full lines) radii versus time. 
4 Momentum Dependent Interaction 

It is well known that the nuclear interaction is momentum dependent. Many 
expression have been proposed, here we use one introduced by C. Gale et al.[6]: 



The values of the parameters are given in [6] . Notice, that also in this case the 
surface term was adjusted to reproduce the binding energy of ^^O. It turns 
out that the value of is half of the canonical one. 

In Figs. 4 and 5 we plot the BE and radii of the nuclei as in Figs. 2 and 3. Also 
in this case convergency is rapidly obtained and the BE are in good agreement 
with the Weizsacker parametrization [4] (dashed lines). The behavior of the 
radii is as before with a neutron halo appearing in the ^^C case. Also notice 
the different oscillations for neutron and protons for this exotic nucleus. 



5 Conclusions 

In this paper we have presented a novel approach to get fermion ground states 
within a Vlasov approach. This approach is the semiclassical counterpart of the 



t/(p(r),p) = ^(p) + -/rfV 



(9) 



1 + [|P-P'|/A]' 
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Fig. 4. Same as in figure 2 but for the momentum dependent interaction. 
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Fig. 5. Same as in figure 3 but for the momentum dependent interaction. 

HF method. It has some advantages compared to HF in a simpler numerical 
solution and a much easier extension to collision problems. On the other hand 
some quantum effects are not included. This is especially important when 
the mean field is not a smooth function. As in all mean field calculations 
symmetries are preserved and because of that, with the use of the mean field 
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parametrization of this paper, nuclei are spherical symmetrical. Correlations 
might break symmetries but these are not contained here. In collision problems 
one can easily include a two body collision term. The approach might be suited 
for calculating low energy phenomena such as fusion, deep inelastic, particle 
production and collective flow at energies up to few hundreds MeV/nucleon. 
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